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FOREWORD 

Robert  Curtis,  Edward  Bensley,  and  Ferrell  Sandy  all  contributed 
important  ideas  during  the  planning  and  programming  of  the  MUSE  Program. 
James  Valentine  and  Augustine  Kish  designed  and  built  the  digital- to- analog 
converter.  Professor  K.  N.  Stevens  of  MIT  provided  very  helpful  criticism 
and  the  use  of  his  sound  spectrograph.  Drs.  F.  S.  Cooper  and  A.  M.  Liberman 
of  Haskins  Laboratory  lent  encouragement  and  valuable  critiques.  The  author 
is  much  in  debt  to  these  and  many  other  individuals. 
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ABSTRACT 

The  MUSE  system,  an  IBM  7090  computer  program  and  associated 
conversion  equipment,  has  been  designed  for  use  as  a  sound  synthesizer. 
Concise  descriptions  of  complex  sounds  including  human  speech  are  converted 
by  the  MUSE  system  into  sound  pressure  waveforms.  The  inputs  to  the  MUSE 
system  are  specifications  of  the  changing  resonance  frequencies  of  multiple 
acoustic  filter  networks  and  of  the  changing  frequencies  and  amplitudes  of  the 
sources  of  acoustic  energy  that  excite  those  networks.  The  output  of  the  MUSE 
system  is  a  sampled  waveform  calculated  for  each  resonance  by  the  solution  of 
a  second-order  difference  equation.  The  results  are  summed  over  a  single 
system  of  resonances  and  then  the  resonance  systems  are  also  added  together. 
The  resulting  string  of  sampled  waveform  ordinates  is  written  in  digital  form 
on  magnetic  tape.  Conversion  to  a  voltage  waveform  is  accomplished  by  use 
of  the  standard  IBM  729  IV  tape  transport  unit  and  a  simple  digital- to- analog 
converter.  Although  the  quality  of  the  sound  is  somewhat  degraded  by  tape 
wow  and  flutter,  acceptable  and  highly  intelligible  speech  has  been  synthesized. 
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SECTION  I 


INTRODUCTION 

A  phonetician,  wishing  to  test  his  model  of  speech  production,  often 
describes  some  speech  utterance  in  terms  dictated  by  his  model  and  then  uses 
this  description  to  control  a  speech  synthesizer.  Provided  the  speech 
synthesizer  is  a  good  one,  the  quality  of  the  resulting  speech  is  a  good 
indication  of  the  sufficiency  of  his  model.  The  model  may  be  too  complicated 
or  may  not  be  limited  by  physiological  constraints,  but  if  high  quality  speech 
can  be  synthesized  according  to  the  model,  it  commands  attention  as  a  possible 
means  of  gaining  insight  into  that  most  complicated  and  refined  biological 
system  —  the  human  vocal  apparatus. 

The  MUSE  system  is  a  computer  simulation  of  a  class  of  sound 
synthesizers  that  have  been  used  with  success  in  testing  theories  of  speech 
production.  In  spite  of  some  loss  of  fidelity,  the  computer  simulation  can  be 
used  in  place  of  this  class  of  sound  synthesizers.  MUSE,  consisting  of  an 
IBM  7090  computer  program  and  simple  digital- to- analog  conversion  equipment, 
translates  concise  descriptions  of  a  large  class  of  complex  sounds,  including 
human  speech,  into  the  corresponding  analog  waveform.  This  signal  can  be 
recorded  for  later  playback  or  it  can  be  used  to  drive  a  loudspeaker  for 
immediate  presentation  of  the  sound. 

Sound  synthesizers  can  also  be  used  to  present  information  in  the  form 
of  spoken  messages  to  the  human  operators  of  a  computer-centered,  real-time 
control  system.  ^  In  such  an  application,  these  messages  would  be  stored  in 
computer  memory  in  concise  digital  form.  The  main  program  for  the  control 
system  would  select  some  appropriate  message  and  send  it  to  a  sound 
synthesizer  that  would  then  expand  the  concise  form  of  the  message  into  the 
corresponding  wide-bandwidth  speech  signal  and  present  it  to  the  operator  in 
real  time.  The  MUSE  system,  by  simulating  the  output  of  various  alternative 
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synthesizers,  could  aid  in  designing  the  simplest  adequate  device  for  each 
application. 

The  MUSE  system  is  itself  a  weak  theory  of  speech  production.  MUSE 
is  a  weak  theory  because  it  leaves  unspecified  the  manner  in  which  the  sounds 
are  produced  and  because,  subject  to  a  limited  bandwidth,  it  can  reproduce  the 
output  of  a  large  class  of  acoustical  devices.  The  limitations  that  make  it  a 
theory  at  all  are  mainly  practical  ones.  MUSE  couldn't  be  used  in  practice  to 
simulate  a  symphony  orchestra  because  a  single  chord  could  require  thousands 
of  statements  in  the  input  language.  General  statements  about  these  limitations 
are  impossible  but  they  will  be  made  explicit  by  the  detailed  description  of  the 
input  language  and  the  operation  of  MUSE  both  of  which  follow. 
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SECTION  II 


EQUIPMENT  NECESSARY  AND  OVER- ALL  SCHEME 
EQUIPMENT  NECESSARY 

The  equipment  needed  for  running  the  MUSE  program  consists  of  a 
standard  IBM  7090  EDPM  with  at  least  two  data  channels  and  at  least  one  729  IV 
tape  transport  unit,  a  special  purpose  digital- to- analog  converter,  a  variable 
bandpass  filter,  and  any  electro-acoustic  transduction  system  or,  preferably,  a 
good  quality  magnetic  tape  recorder. 

OVER- ALL  SCHEME 

In  general,  the  operation  of  the  MUSE  system  consists  of  a  translation 
or  calculation  phase  and  a  subsequent  digital- to- analog  conversion  phase.  In 
the  first  phase,  the  data  cards  containing  the  sound  specifications  are  read  in 
as  needed  and  the  ordinates  of  the  specified  waveform  are  calculated.  These 
numbers  representing  the  instantaneous  pressure  of  the  specified  waveform  at 
successive  small  intervals  of  time  are  stored  in  blocks  on  magnetic  tapes. 

When  calculation  is  completed  for  the  sound  sample  or  utterance  being 
synthesized,  the  waveform  ordinates  are  read  from  the  magnetic  tape  storage 
into  the  computer  for  normalization.  These  data  are  then  written  onto  a  new 
tape  in  which  the  inter- record  gaps  between  the  blocks  of  numbers  are 
eliminated. 

In  the  conversion  phase  the  tape  unit,  an  IBM  729  IV  tape  transport  unit 
on  which  the  final  output  has  been  written,  is  disconnected  from  the  data 
channel  and  connected  to  the  digital-to-analog  conversion  device.  This  device 
then  reads  the  six-bit  (64  levels)  waveform  ordinates  off  the  magnetic  tape, 
sets  them  in  a  buffer  register,  and  converts  them  to  a  voltage  waveform 
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smoothed  by  a  low-pass  filter.  The  resulting  signal  can  be  transduced  by  a 
standard  loudspeaker  system  or  recorded  on  magnetic  tape  as  shown  in  Figure  1. 


MUSE 

STATEMENTS 


DIGITAL  TO 

BUFFER 

ANALOG 

CONVERTER 

MUSE  PROGRAM 


MAGNETIC 

TAPE 


ANALOG  TAPE 
RECORDER 


1.  MUSE  Statements  on  IBM  cards  are  read  into  the  computer. 

2.  The  computed  waveform  ordinates  (A)  are  written  onto 
magnetic  tape. 

3.  The  waveform  ordinates  are  converted  into  a  step  function 
(B)  representing  voltage  levels. 

4.  A  low-pass  filter  smooths  the  step-function  into  an  analog 
waveform  (C). 

5.  The  output  of  the  digital- to- analog  converter  is  a  varying 
voltage  which  drives  a  loudspeaker  or  tape  recorder. 

Figure  1.  Functional  Diagram  of  the  MUSE  System 
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SECTION  in 


SOUND  SPECIFICATION  LANGUAGE  AND  FORMATS 
INTRODUCTION 

Since  the  common  denominator  of  all  sound  synthesis  systems  is  an 
output  waveform,  the  important  differences  between  these  systems  lies  in  the 
method  of  specifying  the  desired  sound  and  in  the  faithfulness  with  which  these 
specifications  are  embodied  in  the  output  waveform.  The  level  of  sophistication 
or,  in  other  terms,  the  degree  of  bandwidth  compression  of  the  system  is 
more  or  less  fixed  by  the  specification  language.  The  MUSE  system's  specifi¬ 
cation  language  consists  of  statements  describing  the  changing  acoustic  spectra 
of  the  desired  sound.  Representing  the  instantaneous  frequency  response  of 
independent  resonating  systems  at  given  points  in  time,  these  statements 
contain  the  resonance  frequencies  and  bandwidths  of  the  several  variable 
resonators  that  make  up  these  systems. 

SOUND  SPECIFICATIONS:  DATA  CARDS 

In  the  description  of  sound  for  MUSE,  resonators  are  grouped  so  that 
several  of  them  can  be  excited  by  the  same  energy  source.  These  groups  are 
called  Spectra.  A  specification  of  the  states  of  the  resonators  in  a  Spectrum 
and  the  excitation  function  for  those  resonators  at  some  point  in  time  is  called 
a  Statement.  Whenever  "Spectrum"  and  "Statement"  are  used  below  in  this 
technical  sense,  they  will  be  capitalized. 

Spectrum  Specification 

The  experimenter  describes  the  instantaneous  frequency  response  of  a 
particular  Spectrum  by  supplying  the  resonance  frequency  and  bandwidth  in 
cycles  per  second  of  each  resonance  contributing  to  that  Spectrum.  More 
explicitly,  the  resonance  frequency  refers  to  the  frequency  of  a  pole  in  a 
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passive  electrical  network  which  is  analogous  to  the  acoustical  system  to  be 
simulated.  The  bandwidth  of  the  resonance  controls  the  real  component  or  the 
attenuation  of  that  particular  pole.  The  value  of  a  resonance  bandwidth  is  the 
difference  between  the  two  frequencies  at  which  the  attenuation  of  the  resonant 
network  is  3  db  greater  than  at  the  pole  frequency  under  the  assumption  that 
the  network  has  only  this  single  resonance. 

Excitation  Source  Specification 

Supplying  the  states  of  resonances,  as  described  above,  can  specify  an 
acoustic  resonating  system  at  a  fixed  point  in  time.  If  enough  points  are  picked, 
or  if  points  of  inflection  are  chosen  and  the  program  is  designed  to  interpolate 
between  them,  a  fairly  complete  description  of  the  continuing  response  of  say 
the  human  vocal  tract,  uttering  speech  can  be  made.  To  excite  the  resonant 
system,  however,  some  provision  must  be  made  for  an  excitation  source. 

In  the  MUSE  language,  energy  is  supplied  to  the  resonant  system  in  the 
form  of  a  train  of  shaped  pulses.  The  pulses,  whose  response  characteristics 
are  fixed  at  the  beginning  of  the  run,  can  excite  the  system  at  periodic  or 
pseudo-random  intervals.  The  mode  of  excitation,  buzz  or  noise,  is  specified 
in  each  Statement.  The  fundamental  frequency  of  the  buzz  source  and  its 
amplitude  are  also  specified.  When  the  noise  option  is  used,  the  repetition 
rate  must  also  be  supplied  since  the  pulses  are  shaped  in  terms  of  decibels 
per  octave  above  the  fundamental  frequency.  The  fundamental  frequency  is 
given  in  cycles  per  second.  The  amplitude  multiplier  is  a  two-digit  number 
specifying  the  relative  amplitude  of  the  source  pulse  (it  is  not  a  logarithmic 
quantity). 

Timing  Specification 

Having  specified  the  Spectra  and  their  excitations  for  single  points  in 
time,  these  points  must  be  fixed  at  particular  times  by  entering  in  each 
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Statement  the  time  in  hundredths  of  a  second  since  the  beginning  of  the  sound 
sample. 

It  has  been  mentioned  above  that  sparse  specification  of  the  variation  of 
a  Spectrum  through  time  would  suffice  if  some  kind  of  interpolation  is  assumed. 
The  translation  phase  of  the  MUSE  program  assumes  that  all  variables  in  a 
Statement  (except  specifications  of  time  and  mode  of  excitation)  change  linearly 
between  Statements.  The  values  of  each  variable  in  successive  Statements  and 
the  time  interval  between  these  Statements  determines  the  rate  of  this  change. 

If  a  variable  has  the  same  value  in  successive  Statements,  that  particular 
variable  remains  constant  throughout  the  interval  between  those  two 
Statements. 

Sorting  the  Specification  Card s 

In  the  process  of  running  the  program,  each  Statement  is  read  in  as 
needed.  Since  interpolation  between  Statements  is  called  for,  it  is  necessary 
that  as  the  calculations  reach  the  "time"  value  of  one  Statement,  the  next 
Statement  in  each  Spectrum  must  be  read  into  the  computer.  The  Spectra  are 
independent  of  each  other  so  ordering  the  Statements  when  specifying  multiple 
Spectra  can  be  fairly  involved. 

Although  it  is  not  necessary  for  running  the  system,  it  is  recommended 
that  a  sorting  field,  not  read  by  the  computer,  be  included  in  the  data  cards. 
This  sorting  field  contains  in  order  a  single-digit  field  that  is  zero  ("0")  for 
the  first  card  in  a  Spectrum  and  a  one  ("1")  for  all  other  cards,  the  time  value 
from  the  previous  statement  in  the  Spectrum,  and  the  Spectrum  number. 

Data  Card  Formats 

The  format  of  the  data  cards  whose  fields  have  been  described  above 
are  as  given  in  Table  I. 
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TABLE  I 


FORMAT  OF  A  SOUND  SPECIFICATION  DATA  CARD 
Field  or  Variable  Name 

Mode  of  excitation  and  card  type  identification:  M0"  *  periodic 
excitation;  "1"  =  noise  excitation 

Number  of  Spectra  ("1"  through  "9") 

Time  (in the  form  XXXX.  XX  sec) 

Amplitude  Multiplier  (XX  arbitrary  units) 

Fundamental  Frequency  (XXXXX  cps) 

Resonance  Specifications  for  Resonance  Number  i  (i  =  1,  2  . . . ,  8),  where 


bandwidth  =  b^  (XXX  cps),  frequency  =  f.  (XXXX  cps): 


Sorting  Field 

"0M  for  first  card  in  Spectrum,  "1"  for  all  others 
Time  from  Previous  Card  in  this  Spectrum 
Spectrum  Number 


Column 


1 

2 


3-8 


9-10 

11-15 


16-18 
19-22 
23-25 
26-29 
30-32 
33-36 
37-39 
40-43 
44-46 
47-50 
51-5  3 
54-57 
58-60 
61-64 
65-67 
68-71 


73 

74-78 

79 


$ 
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Any  resonances  that  are  not  used  can  be  left  blank. 


At  this  time,  a  MUSE  Statement  is  presented  in  its  entirety  on  a  single 
card.  Possible  future  expansions  of  the  input  language  to  multiple  cards  justify 
the  use  of  "Statement"  as  a  special  term. 

CONTROL  CARDS 

There  are  three  control  cards  that  must  be  used  in  the  operation  of  the 
MUSE  program.  These  are  not  acoustic  data  but  serve  to  set  up  parameters 
and  options  for  the  processing  of  the  sound  specifications. 

Parameter  Control  Card 

The  first  control  card  is  identified  by  a  "2"  punched  in  Column  1.  This 
card  specifies  for  the  entire  set  of  data  which  follows  it,  the  number  of  spectra, 
the  number  of  resonances  in  each  of  these  spectra,  the  response  characteristics 
of  the  source  function,  the  sampling  period  of  the  output  waveform  and  a  field 
controlling  program  options. 

The  number  of  spectra  depends  upon  the  number  of  independently  excited 
resonant  systems  desired  (for  speech  synthesis  the  question  of  number  of 
Spectra  is  discussed  under  "Selection  of  Spectrum  Parameters").  In  the 
present  version  of  the  program,  all  spectra  have  the  same  number  of 
resonances,  hence,  only  one  entry  is  necessary  in  the  ”2"  control  card.  The 
response  characteristics  of  the  source  functions  in  all  Spectra  are  the  same 
and  are  specified  in  terms  of  the  slope  of  their  frequency  response  in  db  per 
octave  above  the  fundamental.  The  sampling  period  is  usually  fixed  at  44  x  10 
seconds,  the  program  option  field  at  "00.  " 

The  format  of  the  first  control  card  is  as  given  in  Table  II. 
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Table  II 


\ 

Format  of  MUSE  Control  Card,  Type  "2" 


Field  or  Variable  Name  Column  ' 

Identification  Field  (contains  a  ”2")  1 

Number  of  Spectra  ("1"  to  "9")  2 

0 

Sampling  Period  (in  secs  x  10  ,  usually  "000044"  3-7 

Number  of  Resonances  ("01  to  "08")  8-9 

Source  Characteristic  (in  db  per  octave  x  10  for 
instance,  "060"  is  6db  per  octave)  10-12 

Program  Option  (usually  "00")  13-14 

* 

End  Cards 


There  are  two  kinds  of  control  cards  which  signal  the  end  of  a  set  of 
data.  Any  card  with  a  "4"  in  Column  1  signals  the  program  that  there  are  no 
more  Statements  in  this  sound  segment  and  that  the  run  is  over.  A  "5"  in 
Column  1  indicates  the  end  of  this  sound  segment  and,  in  addition,  sets  up  the 
program  to  accept  the  specifications  for  another  sound  segment  preceded  by 
its  initial  control  card.  When  either  of  these  end  cards  is  read,  the  program 
begins  processing  the  computed  samples  for  the  conversion  phase.  About  one 
second  of  silence  is  automatically  included  on  the  end  of  each  sound  segment  to 
erase  a  section  of  the  output  digital  tape.  The  silent  interval  prevents  spurious 
bits  immediately  following  the  waveform  samples  on  the  output  tape  from 
ruining  the  synthesized  sound  segment  during  the  conversion  process. 
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SECTION  IV 


SPEECH  SYNTHESIS  WITH  THE  MUSE  SYSTEM:  AN  EXAMPLE 
INTRODUCTION 

Although  implicit  in  the  description  given  in  Section  HI,  the  process  of 
using  MUSE  for  synthesizing  sound  is  best  clarified  by  presenting  an  example. 

A  complex  example,  the  speech  utterance  "All's  well  that  ends  well, "  has  been 
chosen  so  that  all  features  of  MUSE  can  be  demonstrated. 

This  example  represents  approximately  2  seconds  of  speech.  It  took 
approximately  200  seconds  of  computer  time  to  synthesize  this  example. 

The  first  steps  in  the  synthesis  procedure  are  common  to  all  so-called 

[2,3,4] 

"terminal  analog"  synthesizers.1  ’  J  First  the  utterance  is  recorded  and 

[5] 

sound  spectrograms  are  made.  Amplitude  sections,  displays  of  the 
Spectrum  of  the  speech  during  some  selected  short  interval  of  time,  can  be 
made  to  clear  up  ambiguous  areas  on  the  time-frequency-amplitude  graph.  An 
additional  useful  datum  is  an  oscillograph  of  the  utterance  with  the  time  scale 
such  that  the  period  of  the  fundamental  and  the  over- all  amplitude  can  be 
measured  throughout  the  utterance.  These  data  for  the  sample  "All's  well  that 
ends  well, "  from  which  an  analysis  of  the  utterance  can  be  made,  are  presented 
in  Figure  2. 

Selection  of  Spectrum  Parameters 

The  next  step  in  the  analysis  procedure  is  to  decide  how  many  Spectra 
will  be  necessary  to  specify  the  sound  sample  to  be  synthesized.  Perhaps  the 
most  important  factor  in  making  this  decision  for  synthesizing  speech  is  a 
heuristic  one;  i.  e. ,  how  many  Spectra  will  best  fit  the  experimenter's  view  of 
the  speech  process?  For  the  example  presented  here,  three  spectra  are  used. 
Each  corresponds  to  a  different  mode  of  operation  of  the  vocal  apparatus.  The 
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Figure  2a.  Sonogram  of  the  Utterance 


Figure  2b.  Oscillogram  of  the  Utterance 


Figure  2 .  Spectrum  of  the  Utterance ,  nAI  I 's  we  1 1  that  ends  wel  I ,  "  Spoken  by  a 
Native  American. 


* 
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first  spectrum  corresponds  to  the  mouth  and  throat  excited  by  periodic 
laryngeal  pulses.  In  other  words,  it  is  used  for  all  voiced  sounds.  The  second 
spectrum  corresponds  to  any  excitation  of  the  vocal  apparatus  by  noisy  sources. 
Fricatives,  aspirates,  and  some  stops  call  for  the  use  of  this  spectrum.  The 
third  Spectrum  represents  the  contribution  of  the  nasal  cavities  to  nasalized 
vowels  and  consonants. 

A  similar  division  of  the  speech  process  into  three  more  or  less  inde- 
pendent  spectra  is  implied  in  the  synthesizers  built  by  Fant  and  K.  N.  Stevens. 

A  more  concise  description  of  speech  utterances  can  be  accomplished  using 
only  a  single  Spectrum  with  a  corresponding  degradation  in  the  quality  of  the 
synthetic  speech. 

The  other  important  decision  to  be  made  at  this  point  in  the  analysis 
involves  the  number  of  resonances  in  the  Spectra.  Three  resonances  are  a 
minimum  number  for  intelligible  speech.  Higher  frequency  resonances 
contribute  to  the  realism  of  the  speech  and  may  be  important  in  realizing  a 
particular  speaker's  voice.  Since  calculation  of  these  extra  resonances  takes 
computer  time,  some  restraint  is  to  be  exercised.  In  the  present  example, 
four  resonances  per  Spectrum  are  specified.  The  resulting  speech  is  highly 
intelligible  but  the  speakers  are  generally  not  identifiable. 

The  slope  of  the  power  spectrum  of  the  excitation  source  must  be 

I  rl 

decided  upon  at  this  point  also.  Certain  theoretical  considerations  by  Fant1 
give  6  db  per  octave  attenuation  as  a  good  empirical  approximation  to  the 
over-all  spectrum  slope. 


* 


These  synthesizers  are  discussed  in  References  [3] 


and  [4]  . 
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The  values  chosen  for  the  number  of  Spectra,  the  number  of  resonances 
per  Spectrum,  and  the  slope  of  the  source  spectrum  are  entered  into  the  "2" 
type  control  card  as  specified  under  "Parameter  Control  Card,"  page  9. 

SEGMENTATION 

Since  MUSE  interpolates  between  parameter  values  on  successive 
statements,  it  is  necessary  to  select  times  for  the  Statements  between  which 
the  Spectrum  parameters  change  only  linearly.  This  is  done  by  close  exam¬ 
ination  of  the  spectrograms  and  oscillograms  of  the  speech  utterance.  The 
segmentation  is  entirely  acoustic  and  has  little  to  do  with  phonetic  divisions. 
An  easy  way  to  place  the  Statements  and  find  parameter  values  is  to  plot  the 
resonance  frequencies  for  each  Spectrum  as  a  function  of  time.  Figure  3 
contains  these  plots  for  the  example  being  considered  here.  The  resonance 
values  can  be  left  constant  or  can  be  changed  linearly  when  the  amplitude 
multiplier  is  zero.  (For  example,  see  times  0.  00  to  0.  06,  1.  98  to  0.  20  in 
Figure  3. ) 

It  can  be  seen  from  Figure  3  that  the  Spectra  are  time-independent; 
that  is,  the  time  segments  in  difference  Spectra  do  not  necessarily  coincide. 
The  values  of  the  various  parameters  can  be  read  from  these  graphs  at  the 
points  corresponding  in  time  to  the  Statements.  The  parameter  values  are 
then  entered  into  a  coding  form  and  punched  onto  IBM  cards  in  the  format 
described  under  "Data  Card  Formats,"  page  7.  After  sorting  the  cards  on  the 
sorting  field,  the  specification  data  are  only  ready  for  calculation.  A  printout 
data  for  these  cards  is  given  in  Table  III  for  the  sample  utterance. 
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??0  2  50  240 


Figure  3.  Tracings  of  the  Sonograms  (Figure  2a)  from  which  the  Input  Data  to  MUSE 
are  Copied.  (The  Vocal  Spectrum  Number  is  "1,"  the  Noise  Spectrum 
is  "2,"  and  the  Nasal  Spectrum  is  "3.") 
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EVALUATION 


The  resulting  synthetic  speech  is  highly  intelligible  if  somewhat 
artificial  sounding.  Recordings  have  been  rather  widely  demonstrated,  with 
almost  total  comprehension  of  the  utterances  reported  by  the  audience.  A 
spectrogram  and  oscillogram  of  the  original  and  synthesized  versions  of  the 
utterance  "All's  well  that  ends  well"  are  given  in  Figure  4. 


% 
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Figure  4a.  Sonograms  of  the  Original  (Top)  and  Synthesized  (Bottom)  Versions 
of  the  Utterance 


Figure  4b.  Oscillograms  of  the  Original  (Top)  and  Synthesized  (Bottom)  Versions 
of  the  Utterance 


Figure  4.  Comparison  of  Sonograms  and  Oscillations  of  Original  and  Synthesized 
Utterance ,  "A  1 1 's  we  1 1  that  ends  we  1 1 .  " 
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SECTION  V 


FLOW  OF  THE  PROGRAM 


INPUT  PROCESSING 

In  order  to  avoid  a  "hardware"  limitation  on  the  number  of  data  cards  or 
number  of  spectra,  the  data  cards  are  read  into  the  computer  as  needed  by  an 
input  subroutine.  When  a  new  card  is  read,  spectrum  parameters  are  set  up 
and  an  interpolation  increment  for  each  Spectrum  parameter  is  calculated. 
These  increments  are  added  after  each  period  of  the  excitation  pulse.  Although 
"continuous"  interpolation  (i.  e. ,  after  each  waveform  ordinate)  would  more 
closely  approximate  the  smooth  changes  in  the  vocal  tract,  substantial  savings 
in  computation  time  accrue  if  the  interpolation  is  lumped  in  coordination  with 
the  source  periods. 

CALCULATION  PHASE 

The  two  (momentarily)  constant  coefficients  of  the  following  second-order 
difference  equation  are  the  recipients  of  the  interpolation  increments: 


where 


and  where 

b.. 

U 

f.. 

y 

At 


A 

B 


=  Axa  a  .  +  Bxa  „  ,  +  S 
t-At  t-2At 


I”  -2tt  b..At 

=  [2e  1J  cos(2?rf„At) 


bandwidth  of  the  i  resonance  in  the  j  spectrum, 
frequency  of  that  resonance, 
sampling  period, 
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=  amplitude  at  time,  t,  and 
S  =  relative  amplitude  of  the  source  function. 

Evaluation  of  this  difference  equation  for  successive  sampling  periods 
results  in  a  series  of  ordinates  of  the  waveform  of  a  single  resonance  excited 
by  a  pulse  of  amplitude,  S,  once  per  period  of  the  source  function.  At  each 
sampling  point,  t,  the  ordinates  of  each  resonance  are  summed,  and  these  sums 
are  added  over  all  Spectra.  The  resulting  over-all  ordinate  is 

l  K 

vIE  v  • 

j  =  l  i=l 

where 

i  =  resonance  index, 

j  =  spectra  index,  and 

i  and  K  =  upper  limits  of  the  indices  (the  number  of 
spectra  and  resonances  given  in  the  "2" 
control  card). 

As  the  calculation  proceeds,  the  waveform  ordinates  are  stored 
temporarily  in  blocks  of  12,  000  on  a  magnetic  tape. 

OUTPUT  PROCESSING 

When  one  of  the  final  control  cards  (with  either  a  "4"  or  "5"  in  Column  1) 
is  encountered,  the  final  output  routine  is  begun.  The  blocks  of  waveform 
ordinates  stored  temporarily  on  tapes  are  read  into  the  computer,  scaled, 
packed  and  converted  to  a  steady  stream  of  six-bit  numbers  written  as  a  single 
record  in  low  density  on  the  output  tape.  If  the  final  control  card  is  a  "4, "  the 
program  is  finished  and  control  is  returned  to  the  monitor  system.  If  the  "5" 
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control  card  terminates  the  utterance,  the  program  clears  storage  space  and 
begins  reading  in  the  control  and  data  cards  for  the  next  utterance. 


v 


ft 
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SECTION  VI 


RUNNING  THE  PROGRAM  AND  CONVERTING  THE  OUTPUT  DATA 
INTRODUCTION 

This  section  describes  the  preparation  of  the  input  deck  for  computation, 
the  running  of  the  program  itself,  and  the  operation  of  the  digital- to- analog 
converter. 

RUNNING  THE  PROGRAM 

The  data  cards,  alter  sorting  on  the  field  described  under  "Sorting  the 
Specification  Cards"  on  page  7,  are  preceded  by  the  "2"  control  card,  followed 
by  the  "4"  or  "5"  control  card  and  inserted  behind  the  "*DATA"  card  following 
the  MUSE  program  binary  deck. 

Tapes  needed  in  running  the  program  are  B5  and  B6  and  A10.  The 
B-channel  tapes  provide  temporary  storage  while  the  final  output  is  written  on  an 
A10.  Under  ordinary  operation  there  are  no  stops  in  the  program.  Any  on-line 
printouts  are  used  by  an  observer  only  to  keep  track  of  the  program’s  operation. 
They  do  not  require  action  from  the  computer  operators. 

OUTPUT  CONVERSION  PROCESS 

While  the  program's  operation  is  quite  routine,  the  output  conversion 
process  definitely  is  not  (see  Figure  1  for  an  over-all  schematic  of  the 
conversion  system).  It  is  advisable  to  make  the  connections  between  the 
converter  and  the  tape  drive  under  the  supervision  of  IBM  customer  engineers. 
Computer  time  is  conserved  if  the  tape  drive  containing  the  digital  output  tape 
is  disconnected  from  the  computer  during  a  routine  halt  between  two  runs.  An 
extra  "terminator"  for  the  detached  tape  drive  must  be  attached  to  the 
converter.  The  bandpass  filter  on  the  output  of  the  converter  should  be  set  at 
a  nominal  low-pass  cutoff  of  about  10,  000  cps. 
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The  tape  recorder  should  be  started  before  the  tape  drive.  A  switch  on 
the  converter  supplies  the  signals  necessary  for  starting  the  digital  tape  drive. 
The  conversion  process  then  continues  in  real  time.  When  conversion  is 
finished,  the  digital  tape  can  be  stored  for  later  use  and  the  tape  drive  can  be 
re- attached  to  the  computer. 

DIGITAL-TO- ANALOG  CONVERSION  EQUIPMENT 

The  digital-to-analog  converter  that  was  used  in  the  MUSE  system  is 
only  one  of  several  circuits  that  would  serve  as  a  conversion  device. 

An  IBM  cable  attachment  plug  is  part  of  the  converter  used.  The  input 
to  the  read  amplifiers  comes  from  this  plug.  The  equipment  is  most 
conveniently  mounted  on  an  ordinary  rack  with  space  provided  for  the  bandpass 
filter  used  to  smooth  the  output  of  the  converter. 
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APPENDIX  A  PROGRAM  LISTING 


Following  is  the  FORTRAN  program  listing  of  the  MUSE  Program. 


L  7] 


H.  SLAHSCN  M  U  S  E  2  ,  A  SPUr  C  SYNTHES  1 /  FR 


A/10/63 


CCMMCN  FRrc.FANC,  SREO,  SANCtCCOA,  CCOO  ,VAl,  TOT  ,  STORE  t  PUT  I  N,  FMI  S, 
xs M  IS t r* Ic  tFCC A, FCCH, FORCE, PER  IOC  ,  I  IN 

X,  I$W2,  J,  SMALL,  I SpFC »  LCOUNT,  T I  ME , 1 R E SON ,  l,  IMS2,  NAR, 
XCRFAT,  SAYPFR 
X  ,  STCR 
X  ,  CL  T  L 

X.PLLSE.OP'Jl  SF.TP 
X , PERTY 
X,  OPTION 

CIMFNSICN  F  R  EC (  8,  10),  F  A  NC (  B,  10),  SRFCI  B,  10),  SAND( 

X  1C),  CCPA(  6,  1C),  DCOBI  ft,  10),  V  A  L  (  2,  8,  10),  TOT (  10), 

X  STCRFI  ICO),  P'jTIM  21)  ,  FMIS(  4,  10)  ,  SMIS(  A,  JO  ), 
XCMS(?tlO),  FCCA(fl,lC),  FCCH  (8*10),  F0RCF(8,10)  , 

X  PFRinC  (  1C  )  ,  I  I N ( 2  ) 

X , STOft ( 12000  ) 

X , CL  T L ( 2CCC ) 

X , PLLSF ( fl,  1C  )  ,C PULSE! 8, 1C) 

X , PFRTY (  10  ) 

FCL'IVALE\CE  (  STCRE,  STOR  ) 


MU2C0300 

MU200AOO 


8 ,MU2C0500 
MU200600 
MU2C0700 

MJ2  00900 


HCISFKEEP  AM;  R  E  AC  CONTROL  CARD. 

MU201100 

ICC 

RFAC  INPUT  TAPE  5,10,1IN(1),  ISPFC,  SAMPFR,  IRESPN,  OB, 

OPT  I CN 

10 

FCRYAT  ( 2  I  l  , F5.0,  I 2, F3. 1  ,  F2.0  ) 

IF  (  CPTICN  -  2.  )  25,15,15 

15 

PRINT  II 

1  l 

FORMAT  (  1FI  ) 

PRINT  12  .  I  IM  1  )  ,  1SPEC, SAHPER, IRESON,  OB,  OPTION 

\ 2 

FORMAT  (IF  12,  13,  E15.7,  13,2615.7) 

25 

VsRlTC  OUTPIT  TAPF  6,  ll 

FRlTE  CUT  PL  T  TAPE  6,  12,  IlN(l),  ISPEC,  SAMPER ,  IRESCN, 
SAMPER  *  SAMPFR  •  10.  **(-6  ) 

06, OPT  ION 

l COUNT  *  0 

MU201300 

SMALL  =  c. 

GPEAT*0. 

RtViINC  8 

MIJ2  0  l  500 

L  *  l 

MU201600 

ISW2  =  2 

NAR  *  I2CCC 

IF  (  CPTICK  -  5.  )  A5,  50,  feC 

MU201700 

45 

IF  f  CPTICN  -  A.  i  60,57,60 

50 

NAR  *  599 

GOTO  60 

57 

N  A  R  =  ICO 

60 

cc  no  j  *  i ,  i  s  p  ec 

PC  112  1  =  1, IRESCN 

MU202000 

V  A  L ( 1 ,  1 , J )  =  0. 

MU202100 

112 

V  A L  (  2  ,  I  f  1)  *  0. 

PERTY(J)  *  0. 

MU202200 

l  10 

TCT(J)  =  0. 

MU202300 

CC  55  1*1,  NAR 
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ooo  r>  o  n  iron  oooooo 


A  .  SL  A  W  SCN 


NUSC2,  A  SQUNC  SYNTHESIZER 


4/10, 


55  STCRl  I  )  =  C. 

TINE  =  0.  MU 2024 00 


END  GF  HOUSEKEEPING.  SET 
UP  INITIAL  CONDITIONS. 


115  CC  130  J=I,  ISPEC 

J- J 

call  SUf* 2  {  INS?  ) 

IJ3  CALL  Si  8  3 

CC  147  JM,  ISPEC 

J  r  J 

CALL  S  U  H  7  (  IVS2  ) 

CALL  SUR4 
14/  CALL  SUP5 

IF  (TPTICN  )  15C,  150,  149 

149  PRINT  4  C  20 

WRITE  CL' T  PL  f  TAP[  6,  4070 
4070  FCRYAT  (  12E  H  E  FOR  E  15C. 

CALL  wCUNP 


PECIN  VAIN  FLCW 


MU202600 


150  J  ^  l 

155  TCTfJ)  O.C 

6C0  DC  620  I  =  l  »  I  RE  SC  N 

TFNP  =  ( FCC A (  I »  J  )  •  VAL ( 1 » I • J  ) )  -  FCCB(I,J) 

X#  VAL (2  ,  I  ,  J  I  ♦  FCRCF (  I  , J  ) 

FCRC F (  I  , J  )  =  C.C 

V  A  L  C  2  •  K  »  J  )  =  V  A  L  (  l .  I  •  J  ) 

V  A  L (  I  ,  I  , J )  =  TEMP 

620  TCTUJ  *  TCT(J)  4  TEMp 

1553  IF(SMKS(  l , J )  -  TIME)  163,  163,1552 
1  552  IF  ( PERT  Y ( J I  -  TIMF  )  160,  169,  157 

157  I  (  ISPEC  -  J  )  159  ,  159,  158 

1 5  fl  J  =  J  4  l 

GC  TC  155 


MU 205300 

MU2 19000 
MU219900 


MU2  2  0200 
MU220300 


MU206 100 
NU206200 


159  TFNP  =  C.C 

DC  820  J*l,  ISPEC 
070  TEM?  ^TCT(J)4  TEMP 
STCR(L)  =  TFNP 
TINF  *  TIME  ,  i. 

L  =  LM 

IF  {  NAR  -  L  )  1591,  150  ,  150 


MU223900 

MU224000 

MIJ2  70500 


1591  CALL  SUP  / 
1595  L  =  l 

GC  TC  150 


MU2O60OO 

MU206900 


PERIOD  OF  THE  FUNDAMENTAL 
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nnonn  r>  n  n  r*  r>  n  nnnnn  o  o  r* 


U.  SLAWSCN  MUSE2,  ^  SOUNC  SYNTHESIZER 


A/  10/61 


IS  OVER.  INCREMENT  DARA^fT- 
ERS  AND  START  A  NEF  PERIOD. 


160 

DC  920 
PULSEt I , J ) 

1  = 

1,  IRFSCN 
PULSEt  I , J) 

♦ 

OPLLSEt I , J ) 

•  PER  I  CO ( J ) 

FCCAII , J) 

- 

FCCAt  I  ,  J) 

♦ 

CCCAt  I#J  ) 

•PERIODt J) 

MU226A00 

920 

FCCHI I.J) 

= 

FCCRt  I, J) 

4- 

DCCRt  I,J  ) 

•PERIHDt J) 

MIJ2  26  500 

F  M  I  S  (  2  ,  J  J 

- 

FMI S( 2,  J) 

4 

OMISt  1, J) 

•  PERI  00  t  J ) 

MU?  2  6600 

161 

F“IS( 3, J) 
call  SUR  5 
GC  TO  157 

FMI St  3,  J) 

4 

cm*  I  S  t  2,  J) 

•PERIOD! J) 

MU?  26  700 
MU207200 
MU207100 

PROCESS  A  NEW  CART 


161  CUL  SOB  1 
1 6 A  C  ML  SCB 2  (  IMS?  ) 

IF  (  IMS?  -  I  )  1 6  7 »  166*  16  5 

167  CALL  SURA 

GC  TC  1552 


Ml)2  0  7  500 
MU207600 
MU 207700 
MU2  08 100 


ENO  CARD  HAS  BEEN  REACHEC. 
START  FINAL  PROCESSING. 


165  CC1651  IX  =  L  *  N  AR 

1653  STCR  ( IX )  =  O.C 

CALL  SLR  7 

IF  (CPTICN  -A.  J  175  ,  1651  ,  1651 
175  re  1751  IX  =  l.NAR 

1751  STCR (  I  X )  =  C.C 

LCCUNT  =  LfCLNT  ♦  6 

CC  1752  I  X  =  1  ,  6 

1752  CALL  RITF(STCR) 

R  E  U INC  P 

R  EW INC  1C 

DC  l R  5  IX  =  1, LCCUNT 
180  CALL  REEC(STCR) 

5002  FORMAT  (  IE  E  19.7,  5E20.7  )  19M9?90 

105  CALL  SUP1C  ( SMALL .GREAT, NAR, S rOR,OLTL ) 

ENDFILC  1C 

1)0  CALL  SURll  19M9290 

1651  PRINT  <.016 

WRITE  CU'iPUT  TAPE  6,  <.016 
AO  1 6  FCRMAT  (12H  ENC  CF  JCR.  ) 

CALL  wcump 

652  I  F ( I  MS  2-2 )  192,192,191 

m  ccto  ioo 

192  CALL  EXIT 


ERRCR  RCL  TINES 
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U.  SLAWSCN  HUSE  2 1  A  SOUNC  SYNTHESIZER 

166  PRINT  A  0  1 7 

ViR  I  T  E  CUT  PL  T  TAPE  6,  A017 

AO  1 7  ECRHAT  (  A  1 E  CATA  CARC  OUT  OF  ORCFR.  STOP  EXECUTION.  I 
CALL  ViCUHP 
CALL  EXIT 

rNC(i,i,o,o,c,o,c,o,o,o,o,o,o,c,o) 
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S  L  H  2  ,  READ  TATA  CARD  INTO  S(J)  REGIONS 


4/10/63 


SUP ROUT  I NE  SGH  2 (  l  R R S  2  ) 

CTMMON  FREQ, F AN C, SRFC,S AND, CCCA, DCOB, VAl, TOT,  STORE, PUT  IN, FM1S. 
XSM  IS  ,r^IS,FCCA  ,FCCB,  FORCE  ,  PER  10  D  ,  UN 

X,  IS  W2  ,  J,  SMALL,  ISPEC,  LCOUNT,  T I  ME  ,  I R  F  SON ,  L,  I  MS  2 ,  NAR, 
XGRFAT,  SAMPFR 
X  ,  STOR 
X  ,  :  l  T  L 

x , PULSE  *C PULSE , CH 
X , PERTV 
x,  CPTICN 

I  I  M  E  N  S  I C  N  FRFCt  8,  10),  F  AND  (  8,  1.0),  SREC(  8,  10),  SANCl 

X  10),  fCOA(  8,  1C),  CCOBt  8,  10),  V  A  L (  2,  8,  10),  TOT (  10), 

X  STORE (  ICC),  PUT  I N (  21)  ,  FMIS(  4,  10)  ,  SM(S(  4,  10  ), 

XrMIS(?,lO),  F  CO  A ( 8 ,  1 0 ) ,  FC0B(8,10),  F0RCE(R,10)  , 

X  PER  ICC  (  1C  )  ,  I  IN ( 2  ) 

X.STCRU20CG) 

X,CLTL (2CCC  ) 

X ,PLLSF ( 8,  10)  , CPU LSF( 8, 10) 

X  ,PFRTY (  10  ) 

FCUIVALFNCE  (  STCRF,  STOR  ) 

200  READ  INPUT  TAPE  6,  2CCC,  IIN(l),  I  I M  2 ) »  ( PUT  1 N ( H ) «  M*l*l<l  ) 

Ir  (  CPTICN-  2.  )  260,  280,  25C 

250  PRINT  2001  ,  IIN(l),  I  I N  (  2  )  ,  (POTIN(M),  M=:l#19) 

2001  FC  R  MAT  (IFC2I2,  F7.,  F4.,  F8.  ,  8(F5.,F6.)) 

260  WRITF  OUTPUT  TAPE  6,  2C01,  IIN(l),  IIM7),  (  PUTIN(M),  M=1,19) 

201  IF  (  I  I M )  )  -  2  )  210,  202,  208 

202  PD  I N T  2L 0 2 

WRITF  OUTPUT  TAPE  6,2002 

2002  FORMAT  (24F  MISPLACEC  CONTROL  CARO.  ) 

CALL  WDl'MP 

CALL  EXIT 

205  IF  (  IIN(l)  -3  )  206,  207,  245 

206  CALL  W  CUM  P 
(ALL  EXIT 

207  PRINT  2C03 

WRITF  OUTPUT  TAPE  6,  2003 

2003  FORMAT  (  50F  CHANGE  IN  NUMBER  OF  RESONANCES.  TO  BE  PROGRAMMED. 
CALL  WCU'WP 

CALL  EXIT 

245  IFIIIMIM)  208,208,240 

208  i  RR  S  2  =  2 
200  RETURN 

C  ENC  CARC. 

210  I  F  (  II N (2)  -  J)  215,  220,  215 

218  IRRS2  *  1 
PRINT  4003 

WRITE  OUTPUT  TAPE  6,  4C03 

4003  FORMAT  (  48F  CATA  CARTS  FROM  WRONG  SPFCTRUM.  STOP  EXECUTION.) 


tall  wcump 

CAUL  DUMP 

220 

S  M  I  S  (  2  *  J  )  = 

PUTINI2) 

225 

S  M 1 S ( 3 ,  J) 

x  1.  /  C 

PUT  I N ( 3 ) 

•  SAMPER 

SMIS(ltJ)  * 

PUT  IN (1)  / 

(ICC.  ♦ 

SAMPER  ) 

8  M  I  $  {  4  ,  J  ) 

V*  4 

*  I  IN  (  1  ) 

MU208800 
MU208900 
MU  209  000 


8 , MU2  09 1 00 
MU209200 
MU209300 

MIJ209500 


MU209700 


) 


MU210400 

MU2 1 0600 
MU2 10700 
MU210R00 


MU2 11400 
MU21 1500 
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SUM2,  «‘Ar  CATA  CAR!  INTC  S(J)  REGIONS. 
N  *  5 


TUC P I  =  ?. 

cr  ?ic  i=i. 

(RESCN 

SREC (  I  , J  ) 

*  PL  T I N ( 

N  ) 

*  TwUPI 

SANCI ( , J 1 

=  PLT(N( 

M  ) 

•TWOPI 

N  =  N  ♦  2 
210  M=M*2 

IRRS2  =  0 
OC  TO  20^ 

C  NLRMAL  RETURN. 

240  ( P  R  S?  =  1 

core  ?C9 

2000  FORMAT  (211  ,  F5.0  ,  F2.0  ,  F6.0, 

t.NC(l,I,O,C,C,0,0,OlO,OtO,C,C,0,O) 


4/10/63 
MU21 1600 
MU21 1700 

MU2U900  f 

MU21 IftOO 

MU212000 

HU212100 

MU212200 

MU212400 


8(  F3.0,F4.0  )  )  MU212500 


SUB3 ,  MCVF  S ( J )  TC  F ( J )  REGIONS. 


4/10/63 


SUPRCUTINF  SLR3 

CCM*CN  RRFCt F ANC» SR FC,SANCtCCOA, ncnB*VALt TOT, STORE » PUTIN, FM(Sf 
X$MIS,CM (S.FCCA, FCCB, FORCE , PER  IOC  » I  I N 

X,  ( Sw2 »  J,  SMALL,  IS^FC,  LCOLNT,  T ( ME  t (RE  SON ,  L,  IMS2,  NAR, 
XGREAT,  SAMPER 
X  f  STOR 
X  ,  CLTL 

X, PULSE, CRUISE, CH 
X, PERTY 
X,  CPIICN 

C  (  MENS  I  GM  F  R  FC  (  8,  10),  F  ANC  (  8,  1C),  SRECt  8,  10),  S  AND  ( 

X  10),  CCOAi  ft,  1C),  CCOfl(  8,  10),  V  AL (  2,  H,  10),  TCT(  10), 

X  STCRM  ICO),  PtTIN(  21)  ,  FMIS(  4,  10)  •  SM(S(  4,  10  ), 

xr'MIS(2,lO),  F  C  0  A  (  ft  ,  10),  FC0H(R,i0),  F0RCE(8,10)  , 

X  PFRICC  (  1C  )  ,  I  IN ( 2) 

X , STOR ( 12900) 

X , OL  T L ( 2 OCC ) 

X .PULSE (8 , 10) t CPULSEt  8,  10) 

X  ,  P  t  R  T  Y  ( 1 0  ) 

ECU  I V  Al C  NCF  (  STORE,  STOR  ) 

300  IS*  M 

CC  310  MM, 4 

310  F  M  I S (  M , J  )  =  SMIS(  M,  J  ) 

CC  320  1=1,  (RESCN 

r a n r (  i , j  )  =  sancc  i,  j  ) 

320  r R  bC (  I,J  )  =  SRECt  I,  J  ) 

RETURN 

CNC ( l, l ,0,0,C,0,C,0,0,0,0,0,0,0,0) 


MU212900 
MU213000 
MU? 13100 


8 , MU2 l 3200 
MU21 3300 
MU213400 

MU200900 


MU 2  l  3700 

MU213900 

MU214000 

MU214100 

MU214200 

MU214300 
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sue  a 


COMPUTE  COEFFICIENT  INCREMENTS 


400 


410 


420 


SUBROUTINE  Sie4 

MY  ENC  IS  MY  BEGINNING.  FIOOLE  WITH  THF  M IDOLE  COEFFICIENTS. 
COMMON  FRFC,FANO»SRFQ»SANC, CCOA ,OCOH , VAL # TOT , STORE , PUT l N , FM I S , 
XS*IS,CMIS,FCCA,FCOB,FORCE,PERIOO  ,  I  IN 

X,  ISW2,  J,  SMALL,  ISPEC,  LCOUNT,  T I M  F , I  RE  SON ,  L»  IMS2,  NAR, 
XGREAT,  SAVPFR 
X  ,  STOR 
X » OL  TL 

X  fPLLSEtCPULSEiCB 
X , °ER  TY 
X,  CPTICN 

C I  ME  NS  I  ON  FREC (  B,  10), 

X  10),  TCOA(  8,  1C),  CCDR ( 

X  STCRE (  100),  PUT  IN (  21)  , 

XCVIS(?,10),  FCOA( 8, 10) ,  FCOBC 8, 10), 

X  PERICC  (  1C  )  ,  I  I N ( 2 ) 

X , STCR ( 12CC0) 

X,CLTL(2000) 

X ,PLLSF (8, 10) ,C PULSE! 8, 10) 

X, PERTY(IO) 

FCUIVALFNCE  (  STCRE,  STOR  ) 

CTIME  -  SMIS(l,J)  -  FMIS(l.J) 

C  =  2.7182818 

CC  410  1=1,  IRFSCN 


FANCI  8,  10) 

8,  10),  V  A  L ( 

F  M I S (  4,  10)  ,  SM I S (  4, 

F  ORC  E (8,10) 


>  SREQ ( 
2 ,  8 , 


8,  10),  SANO  ( 
10),  TOT (  10), 

10  )  . 


FM I S ( 2 , J  )  * 
FMISO.J)  *  SAMPER 


S  1 NF  ( 

)  /  6, 


PULSE (  I, J)  * 

X  *( (  F  REO ( I , J)  •  ( 

X  .  1  6  609ft  ) 

CP  L  LSF (  I  , J  )  *  ( 

X  • ( (  S REO  (  I  , J  ) 

X  .166096  )  ) 

FCC  A ( I , J ) 

X  (  I  ,  J  )  • 

P  C  C  A  (  I  ,  J  ) 

X ( I , J)  • 

FCCeiltJ)  *  E  • •  (-2.0  •  FANO  C  I  ,  J  ) 

CCC8CUJ)  =  (  E  •  •  (-2.0  *  SANO  (  I  ,  J )  < 

X  /  DTIME 

OC  420  1*1,2 

UM  l  S (  l,J) 

RETURN 

ENC  (  l  ,  l,0,G,C,0,C, 0,0,0, 0,0,C,0,0) 


FREOCl , J ) 
2831053))  < 


(  SM  l  S  ( 2  ,  J  )  * 

•  (  SM I S ( 3 , J  )  •  SAMPER 
-  PUL S E (  I , J)  )  /  CTIME 
(2.0  *  E  **  (-FANO(  I  ,  J  ) 
SAMPER  ) 

((2.0  •  E  *  * ( - S ANC  (1,J) 
SAMPER)  )  -  FCO A (  I, J)  )  / 


S 1  N F  (  SREQ  ( l  ,  J  ) 

)  /  6.2831853) )  < 


SAMPER 
(-DB  * 

SAMPER 
(-08  * 


*  SAMPER)) 

*  SAMPER) 
DTIME 


•  COSF  (  FREQ 


COSF  (SRCC 


*  SAMPER) 

SAMPER )  -  FC08 ( I , J ) ) 


(SMIS(IHtJ)  -  FMIS  (1M  ,  J  )  )  /  DTIME 


MIJ2  14700 
MU 2  14900 
MU  2  l  5000 


MU215100 

MU215200 

MU215300 


MU215600 

MIJ215700 


MU  2  1 5  800 
MU215900 
MU216000 
MU216100 
MU216200 
MU216300 
MU216400 
MU2 16500 
MU216600 
MU216700 
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SLB5,  SET  PER  IOH  AND  FORCING  FUNCTION 
SUBROUTINE  SUfi  5 

CCPPCN  FRFC.F ANC.SREO.SAND.DCOA.nCOH , VAL , TOT , STORE , PUT  I N , FM I S , 
X5PIS,CPIS,FCCA,FCCB,FORCE,PERIOC  , I  IN 

X.  ISW2,  J,  SPALL  t  ISPEC,  LCOUNT,  TIME*  IRE  SON  t  L,  IMS2,  NAR , 

XGREAT,  SAPPER 
X  ,  STCR 
X.CLTL 

X, PULSE ,CPULSE,CB 
X , PERT Y 
X,  CPTICN 

C I  PENS  ICN  FRECI  8,  10),  FANCI  8,  10),  SRFOI  6,  10),  SANOI  8.MU217400 

x  io),  rcoAi  a,  ici,  ccoei  a,  ioi,  vali  2,  a,  ioi,  tcti  ioi,  PU217500 

X  STCRE (  ICO),  PUTIN!  21)  ,  FMIS!  A,  10)  ,  SMIS!  A,  10  ),  PU217600 

XCP  IS ( 2 , 10 )  •  FCOAIfi.lO),  FC08 ( 8 , 1 0 ) ,  F0RCE(R,10)  , 

X  PER l CC  (  1C  )  ,  I  IN! 2 ) 

X,  STCR! 1200C) 

X , CUTL ( 20C0 ) 

X,?ULSE(8, 10),CPULSE(8,10) 

X, °ERTY ( 10) 


ECUIVALTNCE  (  STORE,  STOR  I 


IF  (  FP I S ( 4 , J )  -  1.)  500,  5 1 C ,  500 

MU217000 

500 

"ERICC(J)  »  FPISO.J) 

CC  TO  520 

MU2 18100 

510 

CFhRAN  =  RCUN(CUPPY) 

CEWRAN=  (CEWRANM  l./< IC00.*SAMPFR) ) )♦< l./(  5uuu . *SAMPER ) ) 

517 

PERIfC(J)  •  C  EhR  AN 

5  ?0 

PERTY(J)  =  PERTY(J)  ♦  PERICC(J) 

DC  530  1=1,  IRfcSCN 

530 

FORCE  (  I.J)  -=  PULSE!  l,J) 

RETURN 

ENC! 1,1, 0,0,0, 0,0, 0,0,0, 0,0, 0,0,0) 

PU2 1 8  500 

4/10/63 


PU217200 

MU217300 
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SUBROUTINF  7,  WRITE  BUFFER  AREA 


4/10/63 


SUBROUTINE  SUET 

COMMON  FREC,FANC,SRFQ,SAND»CCOA,CCOB , VA L , TOT , STORE , PUT l N , FMI S , 
XSHS,rM$,  FCCAf  FCDB,  FORCE,  PER  100  1 1  IN 

*  X,  l  SW2 •  J,  SMALL,  ISPEC,  LCOUNT,  T 1 MF , I R  FSON ,  L,  [MS2,  NA R f 

XCREAT,  SAMPER 
X  ,  STCR 
X.CUTL 

X, PULSE, CPULSE,CB 
X , PE  R  T Y 
X,  CPTICN 

C l ME  NS l CN  FR  EC (  fl ,  10),  F ANC (  8,  10),  SREO(  8,  10),  SAND!  8 

X  10),  CCOAI  8,  10),  CCOB I  8,  10),  VAL(  2,  8,  10),  TOT(  10), 

X  STORE I  IOC),  PUT  I N (  21)  ,  FMISI  4,  10)  ,  SM I S (  4,  10  ), 

XDM I S  t  2 , 10 ) ,  FC0A(8,  10) ,  FC0B(P,1C),  F0RCE(R,10) 

X  PERIOD  (  1C  )  ,  I  INI  2  ) 

X , STOR I  120CC) 

X ,  l STOR (  12CCC  ) 

X,  CUTL ( 2000  ) 

X  *  PUL  S  E I  8 ,  10)  ,  C  PUL  S  E  I  R  ,  10) 

X  ^PERTY (  10  ) 

E  CU I  VALENCE  I  STORE,  STOR  ) 


X, I STOR 

, IS  rcR ) 

700 

LCCUNT 

*  LCCUNT 

♦  1 

IF  ( 

OPTION  -  4. 

)  705,  707 

705 

CALL 

RITE (STCR) 

GC  TO 

725 

707 

IF  ( 

OPTION  -  6.  ) 

720,  710, 

710 

PRINT 

5001 

PRINT 

5002,  STORE 

720 

WRITE 

OUTPUT  TAPE 

6,  5 CO  l 

0001  FORMAT  (  36F  NAREA  CONSECUTIVE  WORDS  OF  OUTPUT.  ) 

WRITE  OUTPUT  TAPE  6,  5002  ,  STORE 

5002  FORMAT  (  1H  E19.7,  5E20.7  ) 

725  DC  730  L*l,  NAR 

IF  (  GREAT  -STCR(L)  )  740,  750,  750 

740  GREAT  *  STCR(L) 

CO  TO  730 

750  l  F  <  SMALL  -  STCRIL)  )  730,  730,  760 

760  SMALL  *  STCR(L) 

730  CONTINUE 

770  l  F ( SENSE  SWITCH  5  )  775,790 
775  DIFF  =  GREAT  -  SMALL 
CO  780  L  =  l , N  AR 

780  I  STCR I L )  *  III  STURIL)  -  SMALL  )  /  DIFF  )  *  512.)  ♦  256. 

PRINT  7001 

7 C 0 1  r CRM A T I  1 20H  THERE  IS  A  REPEATING  DISPLAY  ON  THE  CRT.  TO  STOP  IT, 
X PUT  SWITCH  5  UP.  AFTER  HALT,  SWb  UP  TO  STOP  TV,  DOWN  TO  SEE  MCRE. 
790  RETURN 

ENCIi, 1,0, 0,0, 0,0, 0,0, 0,0, 0,0, 0,0) 


MU22 l l 00 
MU22  1200 


,  MU2 21300 
M.U2 2  L 4 00 
MU221500 


MU22 1800 


MU222300 

MU222600 
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SUBROUTINE  WCUMP 


4/10/63 


SUBROUTINE  WCUMP 

CO  NON  FREC,FANC,SREO, SANC,  CCOA , DCDB , VAL , TOT , STORE , PUT  I N , FMI  S  , 

XSV IS,  DM  IS, FCCA, FCCB, FORCE , PER  ICC  ,  I  IN 

X,  ISW2,  Jf  SMALL,  ISPFC,  LCOUNT,  TIME  v IRESON,  (,  I  MS2  »  NAR, 

XGREAT*  SAMPER  * 

X  ,  STCR 
X  »  0  L  T  L 

X, PULSE, CPLLSt fTB 
X  f  PERT Y 
X,  C PT1CN 

CIMENSICN  F  R  EC (  8,  10),  FANEM  8,  10),  SRECI  8,  10),  SANO(  8, 

x  10),  rccAi  e,  io),  rcoBi  8,  io),  vau  2,  8,  10),  tch  1 0 ) , 

X  STCRC(  ICO),  PUT  INI  21)  ,  FMIS(  4,  10)  ,  SMISt  4,  10  ), 

XCMIS(2,10) ,  FCC A ( 8 ,  )  0  )  ,  FC  OB (  P ,  10)  ,  F0RCE(R,10)  , 

X  PERICO  I  1C  )  ,  I  IN ( 2  ) 

X  ,  S  TOR ( 120CC) 

X , CUT L ( 200C  ) 

X, PULSE! 8, 10)  ,CP0LSE( 8, 10) 

X , P  E  R  T  Y (  10) 

EQUIVALENCE  (  STORE,  STOR 
WRITE  TUTPUT  TAPE  6,  3C01 

3001  FORMAT  (119E 

X  OPTION  ISPEC 

WRITE  OUTPUT  TAPE  6,  3002, 

3002  FORMAT (  1H  119,  I2C,  F19.7, 

WRITE  OUTPUT  TAPE  6,3021 

3021  FORMAT  (119E 

X  TIME  SMALL 

WRITE  OUTPUT  TAPE  6,3022, 

3022  FORMAT  (1H  I  1 9 ,  1 2 0 , 4F 20 . 7  ) 

WRITE  OUTPUT  TAPE  6,  3025 

3025  FORMAT  (HE  I  INI  X,  J  ) 

WRITE  OUTPUT  TAPE  6,  3026, 

3026  FORMAT  (  1H  119,  120  ) 

WRITE  OUTPUT  TAPE  6,  3013 

3013  FORMAT  (  11E  PUT  IN  t  X 1  ) 

WRITE  OUTPUT  TAPE  6,  3C04, 

3004  FORMAT  (  IE  F19.7,5F20.7 
WRITE  OUTPUT  TAPE  6,  3003 

3003  FORMAT  (  HE  F  RE  C  (  I,  J)  ) 


WR  ITE 

OUTPUT 

TAPE  6,  3  C  04 , 

FREC 

WRITE 

OUTPUT 

TAPE  6,  3005 

3005 

FORMAT 

(  HE 

FANCUtJ)  ) 

WRITE 

OUTPUT 

TAPE  6,  3004, 

FANC 

V.’R  ITE 

OUTPUT 

TAPE  6,  3C06 

3006 

FORMAT 

(  1  IE 

SR  EC (  I , J  )  ) 

WRITE 

OUTPUT 

TAPE  6,  3004, 

SREO 

WRITE 

OUTPUT 

TAPE  6,  3007 

3007 

FORMAT 

(  HE 

SANC(I,J)  ) 

WR  ITE 

OUTPUT 

TAPE  6,  3004, 

SANC 

WRITE 

OUTPUT 

TAPE  6,  3014 

3014 

FORMAT 

(  HE 

F  M  I  S  (  I  ,  J  )  ) 

WRITE 

OUTPUT 

TAPE  6,  3C04, 

FMIS 

WRITE 

OUTPUT 

TAPE  6,  3015 

3015 

FORMAT 

(  1  IE 

SMISt I,J)  ) 

WRITE 

OUTPUT 

TAPE  6,  3004, 

SMIS 

) 

I  IN 


PUTIN 

) 


) 


J  L 

LCOUNT  IRESON) 

(  J,L, OPTION,  I  SPEC , LCOUNT , IRESON  ) 

3120  ) 


IMS2  NAR 

GREAT  OB 

IIMS2,  NAR,  TIME,  SMALL,  GREAT, DB 
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SLR  ROU TINE  WCUMP 


4/10/63 


WRITE  CUT  PL  T 

TAPE  6, 

3016 

3016 

FCRMAT  (111- 

CM  I S  (  I  r 

J)  ) 

WRITE  CUT  PLT 

TAPE  6, 

3C04 , 

DM  IS 

WRITE  CUTPUT 

TAPE  6, 

3017 

3017 

FCRMAT  (Ilh 

FCC A (  I  , 

J)  ) 

WRITE  CUTPUT 

TAPE  6, 

3004, 

FCOA 

WRITE  CUTPUT 

TAPE  6, 

3018 

3018 

FCRMAT  (11F 

F CC B (  I» 

J  )  ) 

WRITE  OUTPUT 

TAPE  6, 

3004, 

FCOB 

WRITE  CUTPUT 

TAPE  6, 

3008 

3008 

FCRMAT  ( 1  IF 

CCC A (  It 

J  )  ) 

WRITE  CUTPUT 

TAPE  6, 

3004, 

OCOA 

WRITE  OUTPUT 

TAPE  6, 

3009 

3009 

FCRMAT  (111- 

CCCB (  It 

J)  ) 

WRITE  CUTPUT 

TAPE  6, 

3C04, 

DCOB 

WRITE  CUTPUT 

TAPE  6, 

3010 

3010 

FCRMAT  (HE 

VAL( X, I 

,  J  )  ) 

WRITE  CUTPUT 

TAPE  6, 

3004, 

VAL 

WRITE  CUTPUT 

TAPE  6, 

3011 

3011 

FCRMAT  (I1E 

TOT ( J  ) 

) 

WRITE  CUTPUT 

TAPE  6, 

3C04  , 

TOT 

WRITE  OUTPUT 

TAPE  6, 

3012 

3012 

FCRMAT  ( 1  IE 

STORE ( L  )  ) 

WRITE  CUTPUT 

TAPE  6, 

3004, 

STORE 

WRITE  CUTPUT 

TAPE  6 » 

3019 

3019 

FCRMAT  (11E 

FORCE ( I 

,  J) 

) 

WRITE  OUTPUT 

T  An  E  6, 

3C04, 

FORCE 

WRITE  CUTPUT 

TAPE  6, 

3020 

3020 

FCRMAT  ( 1  IE 

PER  I  C 0 ( J )  ) 

WRITE  OUTPUT 

TAPE  6 1 

3004, 

PER  1 00 

WRITE  CUTPUT 

TAPE  6, 

3023 

3023 

FCRMAT  ( 1 5 W 

CPULSEt  I  ,  J) 

WRITE  CUTPUT 

TAPE  6, 

3004, 

CPULSE 

WRITE  CUTPUT 

TAPE  6 1 

3024 

3024 

FCRMAT  (I5E 

PUL  SE (  It J) 

WRITE  CUTPUT 

TAPE  6 1 

3004  , 

PULSE 

WRITE  OLT PLT  TAPE  6*  3027 


3027  FCRMAT  (  15F  PERTY(J) 

WRITE  CUTPUT  TAPE  6,  3004,  PERTY 

keturn 

ENC(1*0,0,0,C»0,0, 0,0,0, 0,0»0»0»0) 
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TM-03913 /OOOO/OO/O/OO 


SUBROUTINE  1C  SCALING  AND  PACKING  ROUTINE 

SCALING  AND  PACKING  SUBROUTINE 


CCCO?  ENTRY 


LINKAGF  CIRECTCR 

nccoo  cccooccccccc 

OOOOl  6  264220 ICC6C 


0000? 

0634 

00 

A 

CCC6C 

SURIO 

SXA 

CC003 

0634 

CO 

2 

CCC61 

S  X  A 

COCOA 

06  34 

CO 

1 

CCC6? 

SXA 

CQ0Q5 

0500 

60 

A 

CCC03 

CLA* 

0CCC6 

C  60  1 

CO 

0 

CCC77 

STO 

CCC07 

0500 

60 

A 

CCCO? 

CLA* 

00010 

060  1 

CO 

0 

0CC66 

STO 

00011 

0500 

60 

A 

CCC01 

CLA* 

0001? 

0601 

00 

0 

CCC65 

STO 

coon 

0500 

CO 

A 

CCCOA 

CLA 

000  1  A 

0734 

CO 

1 

CCCOC 

PAX 

00015 

I  OOOOl 

1 

CCC  16 

T  X  I 

00016 

0  6  3  A 

CO 

l 

CCC37 

SXA 

C  00  1 7 

U  500 

CO 

A 

CCCC5 

CLA 

00020 

0734 

00 

1 

CCCOC 

PAX 

00021 

1  OOOOl 

1 

C  CC  ?  2 

T  X  I 

000?? 

0  6  3  A 

CO 

1 

CCC47 

SXA 

000  ?  3 

0634 

no 

1 

CCC27 

SXA 

000  2  A 

1  74C60 

1 

CCC25 

TX  I 

0CO25 

0  6  3  A 

00 

1 

CCC6A 

SXA 

000  ?6 

0  7  7  A 

CO 

? 

C  3  7  20 

A  XT 

000  ?  7 

0600 

00 

2 

5  C  l  ?  1 

L  ?  4 

STZ 

00030 

2  OOOOl 

2 

0C027 

T  I  X 

oonn 

0500 

00 

0 

CCC  66 

CLA 

00032 

030? 

00 

0 

CCC  65 

FSB 

OCO  3  3 

0601 

00 

C 

CCC67 

STO 

0003  A 

0774 

CO 

l 

CCCOl 

A  XT 

00035 

0774 

00 

? 

C  3  7  20 

A  XT 

C0036 

0  7  7  A 

00 

A 

CCC06 

L?? 

A  XT 

00037 

0500 

00 

l 

7746? 

L  ?  1 

CLA 

00040 

030? 

00 

0 

CCC65 

FSB 

000  A  1 

0241 

00 

0 

CCC67 

F  CP 

000  A  2 

0  ?60 

CO 

C 

CCC7? 

FMP 

00043 

-0  300 

00 

0 

CCC75 

UFA 

C0044 

0760 

00 

c 

CCCU 

FRN 

00045 

-0  3?  0 

00 

c 

CCC  7  4 

ANA 

C  00  46 

0767 

00 

A 

CC  C  44 

ALS 

00047 

-060? 

00 

2 

50121 

L  ?  3 

ORS 

00050 

1  OOOOl 

I 

CCC51 

T  X  I 

0005  1 

1  00006 

A 

CCC52 

T  X I 

00052 

-3  00044 

A 

CCC  37 

T  XL 

00053 

2  OOCOl 

? 

CCC  36 

T  IX 

0005A 

♦0776C0002225 

OCT 

00055 

0766 

00 

0 

C2225 

WTBB 

00056 

-O5A0 

00 

0 

CCC64 

RCHB 

00057 

r  06  1 

CO 

0 

COO  5  7 

TCOB 

00060 

0774 

00 

A 

CCCOC 

XSAVE 

A  XT 

0006  l 

0774 

00 

2 

CCCOC 

A  XT 

SUB  1  0 


XSAVt  ,4 
X  S  A  V  E  ♦  I  ,? 

XSAVE  +  ?  ,  1 

3,  A 

NAR 

2,  A 

GREAT 

1.4 
SMALL 

4.4  LOCATION  OF  STOR 

.  1 

*♦1.1.1 

L  2  1  ,  l 

5, A  LOCATION  OF  OUTL 

.  1 

*♦1.1.1 
L  2  3  ,  1 
L  24  ,  1 

*♦1  .  1  .-2000 
101.  I 
2000.2 
OUTL .2 
—  1.2,1 

GREAT  COMPUTE  DIFF 

SMALL 

DIFF 

FND  OF  HOUSEKEEPING 

1,1 

2000,2 

6 . 4 

NARFA ♦ l , l  SCALE  AND  PACK  THIS  RECORD 

SMALL 

DIFF 

SEVSIX 

MAGIC 

MASK  1 
36,  A 
OUTL, 2 
*♦1,1,1 
*♦1,4,6 
L2 1 , 4 , 36 
L  22 , 2  ,  l 

0  7760C002  22  5  SDH  OUTTAP 

OUTTAP 
101 

• 

*  *  ,  4 

*  *  ,  2 


SUBRPUT IKE  10  SCALING  ANO  PACKING  ROUTINE 


00062 

077A  00  l  CCCOC 

A  XT 

**»  1 

00063 

0020  00  A  CCC06 

TRA 

6,  A 

0006  A 

-l  03720  0  A  A  202 

101 

IOCT 

OUTL-1999, ,2000 

0U065 

0  00000  0  cccuo 

SHALL 

PZE 

00066 

o  oonoo  o  cccco 

GREAT 

PZE 

00067 

0  00000  0  cccco 

DIFF 

PZE 

00070 

♦0C00CC0CCCC6 

SIX 

DEC 

6 

00071 

0  00000  0  oocoo 

LOSIX 

PZE 

00072 

♦20677 00 CCCCC 

SEVSIX 

DEC 

63.0 

00073 

♦176ACCCCCCC0 

PFIVE 

OCT 

176A000CQC00 

000  7  A 

♦  COOOC  C  00  007  7 

H  ASK  1 

OCT 

77 

00075 

♦233000000000 

MAGIC 

OCT 

233000000000 

00076 

♦OCOOOIOOOCCC 

ONFCEC 

OCT 

1000000 

00077 

0  OOCOO  0  cccco 

NAR 

PZE 

CCC05 

OUTTAP 

EQU 

5 

7  7  A  6  l 

NARFA 

EQLt 

32561 

5C121 

OUTL 

EQU 

END 

20561 

3 
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3LERCUT  INF  11 


VnR  I  T  t  4  LONG  RECURO 


C  3  7  2  2 

ENTRY 

SUB  l  1 

SUBROUTINE  TO 

WRITE  A  LONG  RECORD 

LINKAGE  CIKfcCTCR 

00000 

OCOOOCOCCCCC 

oooot 

626A2201C 16C 

CCC05 

GUTTAP 

EQU 

5 

CCC05 

F  T  APE 

ECU 

5 

C  0002 

AREA 

BSS 

2000 

C  3  720 

L 

FQU 

2C00 

03722 

06  3 A  00 

2 

C  3  7  5  l 

SUB  1 1 

S  X  A 

EXIT, 2 

03723 

0772  CO 

0 

C  2  2  05 

REWB 

OUTTAP 

0372A 

♦  0  776C  C  0022  2  5 

OCT 

0^7600002225 

SDMB 

OUTTAP 

03725 

-0030  CO 

0 

C  3  7  2  6 

TEFB 

*♦1 

03726 

0762  CO 

0 

C  2  2  2  5 

RTBB 

OUTTAP 

03727 

-05A0  00 

0 

C  3  75  7 

RCHB 

109 

03730 

C 06  l  00 

0 

C  3  7  30 

T  COB 

• 

0373  1 

-0022  CO 

0 

C  3  73 2 

TRCB 

*♦1 

03732 

♦  0776CC00  l; 

2  3  5 

OCT 

07  760000 1205 

SOLA 

F  T  APE 

03733 

0766  00 

0 

C  1225 

WTBA 

F  TAPE 

0373A 

05 AO  CO 

0 

C  3  7  5  7 

RCHA 

109 

03  7  **  5 

0 6A0  CO 

r 

C  3  7  56 

CELAY 

SCHA 

T 

03736 

053A  CO 

•> 

c 

C  3  756 

L  X  A 

T » 2 

03737 

-3  02736 

2 

C  3  73 5 

txl 

•  - 2  »  2  *  ARE  A  +  L-N 

0  3  7  AO 

0762  00 

0 

C  2  2  25 

RTBB 

OUTTAP 

03  7 A  l 

-05A0  CO 

0 

C  3  7  5  7 

RCHB 

109 

0  3  7  A  2 
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